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Abstract 

Here I develop the simplest method in order to evaluate whether or not the Vainshtein mech¬ 
anism can operate for a given set of parameters in a given solution. The method is based on 
the formulation of the mechanism in terms of the Stiickelberg functions given in Int.J.Mod.Phys. 
D24 (2015) 1550022 and arXiv: 1305.0475 [gr-qc]. In such a case, the Vainshtein scale appears as 
an extremal condition of the dynamical metric. If we fix the graviton mass, we can define the 
effective Vainshtein scale. Then for parameters where the Vainshtein scale vanishes or becomes 
smaller than the gravitational radius, the mechanism should be absent. At the other extreme, if 
the Vainshtein scale is finite or infinite, then the mechanism can operate. For consistency, if we 
define the Vainshtein scale as an invariant, then we should expect the effective graviton mass to 
become very large when the Vainshtein mechanism operates. On the other hand, if the mass scale 
tends to zero, then the extra-degrees of freedom are free to propagate. For clarity, here I analyze 
the effective mass behavior for the different type of modes. 

PACS numbers: 04.50.Kd, 04.70.-s, 04.70.Bw 
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I. INTRODUCTION 


The first massive gravity formulation was done at the linear level by Fierz and Pauli [I]. 
This first attempt failed after comparing its predictions with the solar system observations. 
The reason for the failure of the theory at this level is the additional attractive effect due to 
the coupling between the scalar component and the trace of the energy-momentum tensor. 
This coupling remains even at the mass-less limit and it is known as the vDVZ discontinuity 
|2], Then after Vainshtein proposed that it is possible to recover the predictions of GR 
at scales near the source if we consider a non-linear formulation of massive gravity instead 
of the standard linear version 0. Introducing non-linearities however, makes the theory 
pathological due to the appearance of a ghost jl]. In fact, the recovery of GR at scales of 
the solar system was due to the presence of this ghost since the negative kinetic energy term 
reproduced by the ghost exactly cancels the positive kinetic energy contribution coming 
from the scalar component. The details about this physical fact are explained in [5]. Due 
to the presence of a ghost, the non-linear theory was considered pathological. Later de- 
Rham, Gabadadze and Tolley discovered that by tuning in an appropriate way the different 
parameters of the potential, it was then possible to reproduce a ghost-free formulation of 
massive gravity. The essential idea is that the pathological terms of the action, after being 
grouped with higher order terms, are total derivatives and they will never appear in the 
equations of motion. The theory formulated in this way is called dRGT |B]. Inside this 
theory, the recovery of GR at scales of the solar system is due to the non-linearities, relevant 
below the Vainshtein scale which is approximately ry (GM/m 2 ) 1 / 3 . In [7] , the author 
formulated the Vainshtein mechanism in terms of Stuckelberg functions. The Vainshtein 
scale then appeared as an extremal condition of the dynamical metric in unitary gauge. 
This is equivalent to say that the scale is an extremal condition of the massive action. In 
this manuscript, I use this formulation in order to explain for which set of parameters it is 
possible to expect the Vainshtein mechanism to operate and for which set of parameters, 
the mechanism is absent. The theory has generically three free-parameters, namely, the 
graviton mass, and the two free-parameters of the potential. By imposing the Schwarzschild 
de-Sitter background condition, the number of free-parameters is reduced to two. If we fix 
the graviton mass parameter, then we can define an effective Vainshtein scale. For the set 
of parameters where this scale vanishes or becomes smaller than the gravitational radius 
GM , the mechanism must be absent. On the other hand, for the set of parameters where 
the mechanism appears, the effective Vainshtein scale must be finite. If it becomes infinite, 
then the screening mechanism should appear at any scale larger than GM. An alternative 
point of view is to set the Vainshtein scale as an invariant. Then we can define an effective 
graviton mass. At scales where the effective graviton mass almost vanishes, the mechanism 
should be expected to be absent and the massive gravitons can propagate freely everywhere. 
On the other hand, if the effective graviton mass is finite, the mechanism must appear. If the 
effective graviton mass is infinite, then the massive gravitons cannot propagate so far and 
then GR can in principle be recovered. In order to prove the consistency of these previous 
statements, I analyze the effective graviton mass for the different modes separately, namely, 
scalar, vectorial and the purely radial component in a free-falling frame of reference. In 
such a case, the analysis is more transparent and it also takes into account the role of the 
Stuckelberg function when the family of solutions with one free-parameter is analyzed [BJ . 
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II. THE SCHWARZSCHILD DE-SITTER SOLUTION IN DRGT 


The black-hole solutions inside the non-linear formulation were already found in mm- 
The simplest and most generic solution for the spherically symmetric solution, was found in 
[8j. All the solutions satisfying the stationary condition, are defined as: 

ds 2 = G tt dt 2 + G rr S 2 dr 2 + G rt (drdt + dtdr ) + S 2 r 2 dfll, (1) 

where: 

G t t = —f(Sr)(d t T 0 (r, t)) 2 , G rr = -f(Sr)(d r T 0 (r,t)) 2 + 7t4 v 

JW r ) 

G tr = ~f(Sr)d t T 0 (r, t)d r T 0 (r , t), (2) 

and f(Sr ) = 1 — — |A (Sr) 2 , with S being the scale factor which depends on the free- 

parameters of the theory. In this previous solution, all the degrees of freedom are inside the 
dynamical metric. Here we are in unitary gauge and then the fiducial metric is: 

f llI/ dx IJ (lx''' = — dt 2 + dr 2 + r 2 (d9 2 + r 2 sin 2 9). (3) 

The solution ([2]) is equivalent to: 

Q 2 

ds 2 = -f(Sr)dT 0 (r, t ) 2 + + <? 2 rW, (4) 

f{Sr) 

where T 0 (r, t) corresponds to the Stiickelberg function. The non-triviality of T 0 (r, t) contains 
the information of the extra-degrees of freedom nua hd m- In fact, it has been demon¬ 
strated in ra, that To(r, t) represents a preferred time-direction of the theory when the 
extra-degrees of freedom become relevant. The solution Q can be classified in two families. 


1. Family 1: One free-parameter and the Stiickelberg function arbitrary 

For this family of solutions, the Stiickelberg function is also a parameter. In other words, 
there is no rule or constraint telling us how this function should be specified. Then in reality 
we still have two free-parameters, namely, one inside the potential and another represented 
by T 0 (r, t), which has a non-trivial relation with the graviton mass. For this case, the relation 
between the cosmological constant A and the graviton mass, is given by [8j: 

9 

* m 

A = —, 

a 

where a corresponds to the free-parameter of the massive action in agreement with the 
notation given in [8J. The relation between the parameters of the massive action is defined 
in agreement with the condition f3 = a 2 . The result ([5]) is telling us that the graviton 
mass is another free-parameter and its dependence with respect to a must take the form 
m 2 ^ a if we want to keep A as a fundamental constant. Here however, I want to keep 
the Vainshtein scale (not the cosmological constant) as an invariant quantity. In agreement 
with the results obtained in [13], the equivalence principle is not satisfied for this family of 
solutions due to the preferred notion of time in agreement with the function To(r, t). The 
equivalence principle is only satisfied if Tg(r, t) = 0 and this is equivalent to the absence of 
the extra-degrees of freedom for this family of solutions. 
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2. Family 2: Two free-parameters and the Stiickelberg function constrained 


For the family of solutions with two free-parameters and the Stiickelberg function con¬ 
strained, the relation between the cosmological constant and the graviton mass is given by 

m 


A = _ m 2 ( 1 _ l )(2 + a -.|), ( 6 ) 

Depending on the interpretation of the theory, the cosmological constant can be zero if 
we fix the graviton mass as an invariant. If however, what we keep as an invariant is 
the cosmological constant, then the value taken by the graviton mass might change. In 
particular, the relation (3 = (3/4 )a 2 can give us a zero A value or an infinite graviton mass, 
depending on the point of view. Independent of the interpretation used, the Stiickelberg 
function has to satisfy the constraint [8]: 


UoM)) 2 = 


1 - f(Sr) ( S 2 


f(Sr ) 


_ _ rriZ 

f(Sr) 0 


( 7 ) 


A global solution of this previous constraint, is given by the Finkelstein-type form 


rSr 


T 0 (r,t) = St ± 


/(«) 


— 1 1 du. 


( 8 ) 


For the family of solutions with one free-parameter but with the Stiickelberg function ar¬ 
bitrary, in principle there is no clear solution for T 0 (r, t). The arbitrariness of T 0 (r, t) was 
considered as pathological for the perturbative analysis considered in 0. An alternative 
interpretation, was however considered in na The connection of this arbitrariness with 
the Higgs mechanism at the graviton level was analyzed in ra- Note that interestingly, for 
the family of solutions analyzed in this section, the equivalence principle is recovered. In 
fact, from eq. ft it is easy to demonstrate that Tg(r,£) = 0 if f(Sr ) —> 1. Then, for the 
first family of solutions (one free-paramater), we can conclude that the equivalence princi¬ 
ple is in reality hidden. In other words, the principle which is absent from the family of 
solutions with one free-parameter and T 0 (r, t) arbitrary, is recovered when we select as free- 
parameters those inside the massive action [13J. Note that the Stiickelberg function depends 
on the graviton mass through the function f(Sr ) defined previously. Then the arbitrariness 
on T 0 (r,t ) defined for the previous family of solutions, is equivalent to an arbitrariness on 
the graviton mass value. Then for the first family of solutions, the two free-parameters are 
the graviton mass and a. 


III. THE VAINSHTEIN CONDITIONS 

The Vainshtein conditions were derived originally by the author in |7J. They express the 
Vainshtein scale as an extremal condition of the dynamical metric. This result helps us to 
find a relation between the Stiickelberg function T 0 (r,t ) and the function f(Sr) contained 
inside the components of the dynamical metric. The key relations are: 


d r T 0 (r, t) = 0 —>■ r << r v 
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<9 r T 0 (r, t) ^ 0 -» r >> r v 
To(r,t) = 0 -a r = r v , 


( 9 ) 





with ry being the Vainshtein scale. The vanishing condition Tj-'/r. t) = 0, guarantees the 
recovery of GR in the usual Schwarzschild-like coordinates. The same set of relations, could 
also be derived from the most direct condition dU(g, 0 ) = 0 [7J. The result is independent 
of the exact form of the potential U(g, <fi). 


IV. THE VAINSHTEIN CONDITIONS: EXTREMAL CONDITIONS ON THE 
MASSIVE ACTION 


The Vainshtein scale, corresponding to an extremal conditions of the dynamical metric, 
can also be derived directly from the potential U(g,(f )) as has been explained in the previous 
section. The key point is to understand that in unitary gauge, namely, when all the degrees 
of freedom (Nambu-Goldstone bosons) are inside the dynamical metric (eaten up by the 
dynamical metric), then any variation of the potential, is directly related to the variations 
on the dynamical metric. The relevant mathematical condition is [7j: 




d(f>. 


( 10 ) 


The importance of the Vainshtein conditions is that the details of the potential are not rele¬ 
vant at all in order to derive the Vainshtein scale. All the information about the Vainshtein 
scale is contained inside the dynamical metric if it contains all the degrees of freedom. For 
stationary cases, where the metric components are time-independent, then eq. (10) becomes: 


dU(G) = 




( 11 ) 


and this is just equivalent to: 


dg = 0, 


( 12 ) 


if the matrix is non-singular. The condition (12), has to be satisfied by each compo¬ 

nent of the metric containing the information of the extra-degrees of freedom. The conditions 
are valid at the Vainshtein scale. Then whenever we want to calculate the Vainshtein radius 
of the theory, all what we have to do is to translate all the degrees of freedom to the dy¬ 
namical metric and then apply the conditions (|12|) and finally solve algebraically for r = ry 

m 


V. THE PARAMETER-DEPENDENT VAINSHTEIN SCALE: THE CASE WITH 
ONE FREE-PARAMETER WITH A FIXED GRAVITON MASS 

Here I derive the effective Vainshtein scale for the two type of solutions obtained in 
[8j. The first solution to be considered is the one with one free-parameter but with the 
Stuckelberg function arbitrary. In this case, in principle there is no correspondence between 
T 0 (r,t ) and the metric function f(Sr). The dynamical metric in unitary gauge takes the 
explicit form given by eq. ([ 2 ]) but taking into account that T 0 (r,t) is in principle arbitrary. 
Then although the location of the Vainshtein scale in terms of the Stuckelberg function is 
not clear at all, we can assume that the Vainshtein conditions can still be applied to the 
dynamical metric, providing then the information about the location of the scale after which 
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FIG. 1: The effective Vainshtein scale for a fixed graviton mass. Note that the (effective) Vainshtein 
scale becomes asymptotically infinite when a 0 and a —» oo. The Vainshtein scale has a 
minimum for a = 2. The factor has been normalized to be one. This figure also shows 

the qualitative behavior of the Vainshtein scale for one of the branches of the solution with f3 = 0 
and a arbitrary. 


the effects of accelerated expansion have to be taken into account. For the stationary case, 
in agreement with j7jj, we can define the effective Vainshtein scale as: 


f'(Sr) = 0, r = ?v, 

1 /3 GM \ 1/3 _ a + 1 / 3GM \ 1/3 

S \ A J a 2 / 3 \ m 2 J 


(13) 

(14) 


This is an effective value of the Vainshtein radius depending on the free-parameters of the 
theory. In this case however, the possible effects of the Stiickelberg function (T 0 (r. t)) have 
not been taken into account. Later in this manuscript I will analyze the role of T 0 (r, t) but 
inside the definition of effective graviton mass. The analysis done in Sec. flVHl ) , takes into 
account the effective mass for each mode if we keep the Vainshtein scale as an invariant. Fig. 
0 shows the behavior of the scale (14) as a function of a. Note that the Vainshtein scale is 
non-zero for any positive value taken by the parameter a. This means that the Vainshtein 
mechanism will always appear for this solution. When ry —>■ oo, GR should be recovered 
everywhere. However, if we analyze term by term in the function f(Sr), given in terms of 
a by: 


f(Sr) 


2 GM(a + 1) 
ar 


1 a 2 r 2 
3(a + l ) 2 ’ 


(15) 


then it is clear that GR is recovered everywhere for a —>■ oo. However, for a —>■ 0, the 
Newtonian term 2Gj 'V +1 ^ becomes dominant, rescaling the effective Newtonian constant 
and making gravity very strong and attractive. Then the effective Schwarzschild radius 
(gravitational radius) becomes very large. This means that even if the Vainshtein mechanism 
operates everywhere outside the gravitational radius, this part of the solution is unphysical 
from the experimental point of view. However, the main point here is that as far as the 
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FIG. 2: The behavior of the scale factor S' as a function of the free-parameters of the theory. This 
plot corresponds to the negative root square branch in agreement with eq. (16). This picture is 
the mirror image of the positive root square branch of S. 


Vainshtein scale is non-zero, the mechanism is able to operate outside the gravitational 
radius. The case explored here is equivalent to the one with the condition (3 = a 2 . 


A. The case with two free-parameters and the Stiickelberg function constrained 
and a fixed graviton mass 


For this case, there is a direct correspondence between T 0 (r, t) and f(Sr). However, still 
we can define the Vainshtein scale in the same way as in the previous case. It is given by the 
result (14), but with the cosmological constant now defined as in eq. (|6]). The scale factor 


S is given by [E] : 


c _ ct + /3±y/a; 2 — /3 
1,2 “ 1 + 2a + p 

The solution can then be classified in terms of the different values taken by the parameters 
a and /3. Figure shows the behavior of S' as a function of the two free-parameters of the 
theory for the negative root square branch. The same factor but for the positive branch is 
just the mirror image of this figure. In what follows, I will analyze two cases. 



B. Case i). (3 = 0, a arbitrary and a fixed graviton mass 

For this case, there are two possible values for the scale factor S. They are: 


Si = 


2a 

l + 2a 


S 2 = 0. 


(17) 
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The cosmological constant as a function of the graviton mass is in agreement with eq. a 
the following: 


2 2 

Tin mm 

Ai = (Si -1)(25! +aS 1 -a), \ 2 = - (18) 

Note that A 2 diverges if a; ^ 0. If a = 0, then we recover the special case of the previous 
section, namely, with one free-parameter but with the Stiickelberg function arbitrary. This 
case is also divergent for a —> 0. Here we can define the parameter-dependent Vainshtein 
scale as the extremal condition of the dynamical metric in unitary gauge [Zj. The result is: 

2a+ 1 / 3GM\ 1/3 / 3GM \ 1/3 

VV1 ~ (6a 2 ) 1 / 3 y~m 2 ”) ’ ry2 “ “ ^ 

It is evident from the result ( [l7| ) that ry 2 will diverge for any value of a. This means that the 
screening effects are expected to appear everywhere up to the gravitational radius scale in 
this particular situation. However, it is easy to find that the effective gravitational radius is 
also divergent, then the Vainshtein mechanism and any physically relevant solution related 
to this branch is absent. For the case of ryi, the behavior is just the same as in the case 
studied previously when (3 = a 2 . Fig. Q still represents the behavior for this case in the 
qualitative sense. In fact, whenever we fix the parameter (3, independent on its relation with 
a, the qualitative behavior of the effective Vainshtein scale will be the same. 


C. Case ii). Both, (3 and a arbitrary and a fixed graviton mass 


For the two branches of solutions defined in agreement with eq. (16), the parameter- 
dependent Vainshtein scale for this case becomes: 


rv 


3GM{1 + 2a + (3) 3 (3 2 


1/3 


a 


i \J a 2 — /3 T f3 I l m, 2 (—2a 3 i 2 ot 2 \J a 2 — f3 + 3a/3 2 /3-\J a 2 — /3) 


. (20) 


The two branches of solutions can be observed in Figs. ©, @ and Note that the effec¬ 
tive Vainshtein scale has a discontinuity region. In agreement with the plots, the Vainshtein 
mechanism operates whenever we have a finite value for ry. 


VI. THE PARAMETER-DEPENDENT GRAVITON MASS: THE VAINSHTEIN 
RADIUS AS AN INVARIANT 


A. Family 1: One free-parameter and the Stiickelberg function arbitrary 


If we fix as an invariant the Vainshtein scale, then we can construct an effective graviton 
mass. In this situation, we can interpret any of the two free-parameters, namely a or /3 as 
the graviton mass. Here I will analyze the cases exposed in the previous sections but by 
keeping the Vainshtein scale as an invariant and then observing the behavior of m 2 . For the 
case with one free-parameter corresponding to the Vainshtein scale given in eq. (14), the 
effective graviton mass should behave as: 














FIG. 3: The parameter-dependent Vainshtein scale for a fixed graviton mass and arbitrary param¬ 
eter combination. The peaks of the figure represent the largest possible Vainshtein scales. There 
is a region of discontinuity around the peaks representing the absence of the mechanism for the 
relevant parameter combination. The factor ) 1 ^ 3 has been normalized to one. The yellow plot 
corresponds to the positive root square branch in agreement with eq. (16) and the blue graphic 
corresponds to the negative root square. 



FIG. 4: The same Figure (§, but from a different perspective. From this image the plane /3-a can 
be appreciated. 
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FIG. 5: The same Figure Q, but from a different perspective. 



FIG. 6: The parameter-dependent graviton mass as a function of the free-parameter of the theory. 
In this case, the Vainshtein scale is kept as an invariant. This case corresponds to the condition 
a = /3 2 


m 2 ^ 


(a + l) 3 

a 2 


( 21 ) 


Fig. @ shows the behavior of m 2 for the different values taken by the free-parameter 
of the theory. Note that the effective graviton mass goes to infinity near a = 0. In the 
neighborhood of this value, the massive gravitons cannot propagate and then GR can in 
principle be recovered. This is consistent with the ideas expressed in the previous section 
suggesting that for a fixed graviton mass, infinite values of the effective Vainshtein scale 
correspond the range of parameters where GR can in principle be recovered. Note that in 
this case we are not taking into account the role of the Stiickelberg function. This situation 
will be analyzed later in this manuscript. 


B. Family 2: The parameter-dependent graviton mass for f3 = 0, arbitrary a and 
a fixed Vainshtein scale 


If we fix the Vainshtein scale as an invariant, the parameter-dependent graviton masses 
corresponding to the case with f3 = 0 and a arbitrary are in agreement with eqns. (19): 
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FIG. 7: The parameter-dependent graviton mass as a function of the two free-parameters of the 
theory. In this case, the Vainshtein scale is kept as an invariant. The yellow plot corresponds to 
the positive brach of S and the blue color to the negative one. 



FIG. 8: The parameter-dependent graviton mass as a function of the two free-parameters of the 
theory. This plot corresponds to the Fig. 0, but observed from a different perspective. 


m? ^ 


(: 2a+iy 


mi, 


cr 


S 2 a 


( 22 ) 


For m 2 , S? —> 0, which indicates that the parameter-dependent graviton mass diverges in 
that case for any value of a. This again means that the massive gravitons cannot propagate 
for this branch of solutions and then GR could in principle be recovered everywhere. How¬ 
ever, as has been discussed before, the gravitational radius for this solution also diverges, 
making it unphysical. The behavior of m i is exactly the same as in the case indicated in eq. 
(21) and then Fig. is also appropriate for representing this branch of solutions. This is 
expected because whenever we fix the parameter /?, the behavior of the solution with respect 
to a will not change in the qualitative sense. 
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FIG. 9: The parameter-dependent graviton mass as a function of the two free-parameters of the 
theory. This plot corresponds to the Fig. Q, but observed from above. 


C. Family 2: The parameter-dependent graviton mass for arbitrary /? and a pa¬ 
rameters with a fixed Vainshtein scale 

In this case, the behavior of the parameter-dependent graviton mass is more complicated, 
but it can be visualized in Figs. 0, 0 and The regions where the parameter- 

dependent graviton mass takes large values, correspond to regions where in principle the 
massive gravitons can travel shorter distances. This is equivalent to an screening effect. 
Whenever the effective graviton mass is almost zero (but still finite), then the massive 
gravitons are able to travel larger scales and the screening effects are absent. However, 
although the method proposed here is simple, it does not take into account the possibility of 
having different behavior for different modes, namely, scalar, vector and tensor. This issue 
will be analyzed later in this manuscript after taking into account the perturbations of the 
action in a free-falling frame of reference. 


VII. PERTURBATIVE ANALYSIS 

In the previous sections, I worked around the background solutions without any concern 
about the perturbations. From the physical point of view however, if we are worried about 
the propagation of gravitons around some specific background, then we should analyze 
the action at the perturbative level. This provides a cleaner physical analysis and in this 
situation it is possible to understand the role of the Stuckelberg function. The analysis at 
the background level, does not provide the opportunity for understanding the role of this 
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function in the propagation of the graviton field. It has been demonstrated recently by the 
author that T 0 (r, t) can be considered as a preferred time-direction, breaking potentially the 
Lorentz symmetry for a free-falling observer [13J . In this section, I consider the perturbations 
of the action in agreement with a free-falling frame of reference. The method was developed 
in [12|, based in the generic perturbation analysis done in [S]. The purpose here is to find 
the effective mass for each of the modes, namely, scalar, vectorial and radial one. The most 
important mode is the scalar one because at the decoupling limit it generates the extra- 
attractive effect if the Vainshtein mechanism is absent. In general, the different modes of 
the graviton field, have different values of mass. 


1. Family 1: One free-parameter and the Stiickelberg function arbitrary 

This case was considered the most interesting one in [13], when the spatial derivative of 
the Stiickelberg function is non-vanishing. In this case, as has been explained in Sec. |n]), the 
equivalence principle is violated due to the preferred time direction which is a physical effect 
of the extra-degrees of freedom. The effect can also be observed as a spatial-dependence of 
the function T 0 (r, t). Under the assumption of stationary solutions, the potential expanded 
up to second order is given by: 


V=gU(9,<l>) 


1 + l h 


\h a y n + lh 2 


1 + l h 


h 


(1 + ol) 


+ 


/ 3 ,L a 

2T'(r,f)(l + a) 2 


2 T 6cy(l T oi) 


flc\r 


(l + ay 

T'(r,f) 2 (l + a) 3 


cr 


or 


+ ... 


(23) 


The matrix mass is obtained after finding the second derivative of the previous potential. It 
is given explicitly by: 


mF v m a[3 = ( + -fj^fj a0 

2 4 


2 + 6a(l + oQ \ -^-a/3 

a+«)* J v 1 


a 


- (fTS^S *o + f, ap r 0 5\) + Q ) 2 5 (24) 


where corresponds to the inverse of the dynamical metric in a free-falling frame of 
reference as has been defined in [13] . It is defined as: 


ds 2 = S 2 —dt 2 + dr 2 


1 - 


vrr) 

s 


_2^Mdtdr + rW 


(25) 


Without solving for the eigenvalues, we can understand the behavior for the graviton mass 
modes. The matrix (24) can in principle be diagonalized. However, it is not necessary at 
the moment of analyzing the different modes. We can find the value of the graviton mass 
for the different modes from eq. (24) if we download the indices by using the background 

(26) 


metric. The general expression is given by: 


7/3 


Then the effective masses for the different modes are: 
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FIG. 10: The square of the effective mass for the scalar mode. It can be observed that for a 
vanishing value of Tg(r, t) (M in the figure), the scalar component of the graviton mass vanishes. 
Different branches appear for this solution. 



FIG. 11: The same Fig. (10) but observed from behind. Here it is clear that the existence scalar 
mode mass parameter vanishes when Tg(r, t) vanishes as a parameter. Tg(r, t) is represented by M 
in the figure. 


m m = S 4 m 00 + S 2 Tq(t, t) 2 m rr + S 3 Tq(t, t ) (m 0r + m r0 ) , 


(27) 


mo? — t)m°° — STq(r, t) (5 2 — T^{r, t ) 2 ) m rr +S 2 T^{r, t) 2 m r0 — S 2 (S 2 — T^{r, t ) 2 ) 

(28) 


m rr = S 2 Tq(t, t) 2 m 00 + (S 2 - T^r, t) 2 ) 2 m rr - ST* 0 (r, t ) (,S 2 - T'(r, t ) 2 ) (m 0r + m 0r ) . (29) 


The figures (10), (13) and (16) represent the behavior of the masses corresponding to 


different modes. It can be observed that the scalar and radial components of the effective 
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FIG. 12: The same Fig. (10) but observed from the left. From this angle, we can observe the 
scalar component of the parameter-dependent graviton mass as a function of a. From the graphic, 
it is clear that the effective mass for this mode can become very large when a 0. 



FIG. 13: The parameter-dependent square of the graviton mass for the vectorial component. In 
this case, some divergences appear when a —>> 0 independent on the values taken by Tg(r, t). 



FIG. 14: The same Fig. (13) but observed from behind. 
Tg(r, t) can be appreciated. 


Here the dependence with respect to 
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FIG. 15: The same 
can be appreciated. 


Fig. (13) but observed from the left. Here the dependence with respect to a 



FIG. 16: The square of the effective mass for the r — r component. Here we can observe that the 
effective mass for this mode diverges for a 0 independent of the values taken by Tg(r, t) taken 
as a parameter. 



FIG. 17: Fig. (|16j) 
and negative values 


but observed from behind. The plot shows a discontinuity when Tg(r, t) —0 
of a. 
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FIG. 18: Fig. (16) but observed from the left side. Here the dependence with respect to a can be 
appreciated, again the divergence when a —> 0 appears for this case. 


mass, have similar behavior with respect to the parameter a and the function Tg(r, t). From 
the Figs. (12), (15) and (18), it is observed that the effective masses for the different modes 
become very large if the parameter a becomes small, independent on the value taken by 
the arbitrary Stiickelberg function. These results are consistent with those obtained in Fig. 
(JgJ) , where it was defined an effective graviton mass, representative for all the modes. In 
such a case however, the role of T 0 (r, t) was not clear because the analysis only covered the 
behavior of the parameter a. The results are also consistent with the definition of parameter- 
dependent Vainshtein scale and as a consequence with the Fig. 0 - In other words, a large 
effective graviton mass, implies the no propagation of the corresponding massive modes if 
we select as an invariant the Vainshtein scale. On the other hand, if what we select as 
an invariant is the graviton mass modes, then an infinite effective Vainshtein scale is also 
equivalent to the no propagation of massive modes. In the case under study, it is lucky 
that all the modes have a similar behavior for some set of parameters. However, in general 
different modes have different effective masses depending on the selected values for a and 
To(r, t) if we keep the Vainshtein scale invariant. 


2. Family 2: Two free-parameters and the Stiickelberg function well defined. The equivalence 
principle recovered 

As has been mentioned before, for the set of solutions with two free-parameters and 
T 0 (r, t) well defined, the equivalence principle is recovered. The potential expanded up to 
second order in a free-falling frame of reference is given in this case by: 


V=gU(g, cf>) ~ S 4 U(g, fi) back + h (y U(g, cf>) back + S 4 F(a, /?)) - y h a ^ a U(g, fi) back 

+h 2 (jU{gA) back + yi>,/3)) -(30) 

If we repeat the same procedures of the previous case, then we find that the matrix mass is 
given by: 

m^m aP = -y U back (g, </>)fTf} vP + S 4 fj^fj a ^ (hj{g, 4>) back + F(a, /3)^j . (31) 
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FIG. 19: The square of the parameter-dependent graviton mass for the scalar mode and the radial 
mode. The yellow color corresponds to the positive branch of the scale parameter S and the blue 
color corresponds to the negative branch. 


Here we can still use the same expansion given in eq. 
the background metric now is different. In fact, for a 
dynamical metric is simplified to be: 


(26), but taking into account that 
free-falling frame of reference, the 


ds 2 = S 2 ds 2 M , (32) 

which is conformal to Minkowski. Then in this case, the definition of the mass for the 
different modes is trivial and given by: 


m^ = S^ri^rivpm 113 , 

with rg tl/ given by the the standard Minkowski metric, 
following results: 


From eq. 


(33) 


(31), 


we can find the 


1 


{m ) = (m rr ) = --U(g,<j>) bac k + F(a,P), 


(34) 


(m r ) = ~U(g,(j)) back . 


(35) 


The masses for the different modes can be obtained from eq. (33) by using the results 


and (35). Figs. (19) and (23), illustrate the behavior for the masses of the different scalar, 


radial and vectorial modes. Note that the scalar and radial modes have exactly the same 


behavior in agreement with eq. (33). If we analyze the plane a-/3, then we can find similar 
behavior with respect to the case analyzed in Sec. (VI) in Fig. (|9]) as can be perceived the 
graphics (22) and (26). However, as can be expected the behavior of the effective mass for 


the different modes differ. 

The behavior for each of the modes with respect to the parameter a or /?, reveal the 
existence of some peaks as can be seen from Figs. (20) and (21) for the scalar and radial 
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FIG. 20: The same Fig. (19) but observed from the frontal perspective. Here the variation with 
respect to a can be appreciated. 



1000 


0 


1000 


FIG. 21: The same Fig. (19) but observed from the left. Here the variation with respect to (5 can 
be appreciated. 


modes and from Figs. 


(24) and (25) for the vectorial one. For the scalar and vectorial 
The peaks in such a case, can be interpreted as 


modes, the observed peaks are very large, 
the combination of parameters for which the mode becomes heavy and does not propagate 
so far. This is equivalent to a screening mechanism effect. In addition, we can observe that 
the peaks for the vectorial component are not so large. This means that the vectorial modes 
in this case can propagate freely for most of the parameter combinations of the theory. The 
vectorial modes however, normally are decoupled from matter when we analyze the theory 
at the decoupling limit level [6]. 
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FIG. 22: The same Fig. (19) but observed from above. The dependence with respect to the two 
free-parameters is similar with respect to the case analyzed in Sec. (VI) as can be perceived in 
Fig. 



FIG. 23: The square of the mass for the vectorial component. The yellow color plot corresponds 
to the positive branch of the scale parameter S and the blue color corresponds to the negative one. 
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FIG. 24: The same figure (23) but observed from the frontal perspective. Here the variation with 
respect to a can be appreciated. 



FIG. 25: The same figure (23) but observed from the right. Here the variation with respect to /3 
can be appreciated. 



FIG. 26: The same figure (23) but observed from above. Here the dependence with respect to 
both, a and /? can be appreciated. The figure is qualitatively similar to the scalar and radial mode 
case analyzed in Fig. (22). 
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VIII. SUMMARY 


Here I have derived the simplest method in order to predict whether the Vainshtein 
screening mechanism might appear for some specific parameter combinations of the theory. 
The method is based on the Vainshtein conditions derived by the author in [7j. The 
details of the potential are not important for the present analysis. The regions where the 
Vainshtein mechanism operates are not necessarily relevant from the the physical point of 
view. This means that the conditions derived in this manuscript are necessary but not 
sufficient in order to guarantee the recovery of GR. The recovery of GR not only implies 
the predictions at the solar system scale but also the predictions related to the physics 
of black-holes. In ra for example, it was discovered that some parameter combinations 
can reproduce branch point effects when we analyze the periodicity structure of the 
propagators. This manuscript is developed by using two equivalent points of view. The first 
one is worked by fixing the graviton mass and then deriving an effective Vainshtein scale. 
The second point of view is developed by fixing the Vainshtein scale and then deriving 
an effective version of the graviton mass. Here I consider that both points of view are 
physically equivalent and considering one or the other is a matter of choice. Finally in 
order to perform a rigorous analysis, I have studied the behavior of the effective mass for 
each of the modes, namely, scalar, vectorial and radial. In this case I have also maintained 
the effective Vainshtein scale as an invariant. The analysis revealed that the mass for the 
different modes can behave in a different way depending on the parameter combinations 
considered. In other words, although from the qualitative point of view some cases look 
similar, the detailed study revealed that the effective mass for the different modes are not 
the same in general. In particular the scalar and radial modes have the same effective 
mass for the case of two free-parameters and T 0 (r, t) constrained. In such a case, only 
the vectorial component behaves different. However, if we look at the case with one 
free-parameter and T 0 (r, t) arbitrary, then the different modes behave different when we 
look at them in detail. However, even in this situation, the scalar and radial components 
have similar behavior. The difference comes from the behavior with respect to the arbitrary 
Stiickelberg function. In fact, this function defines the direction of the preferred notion of 
time. It would be interesting to extend these ideas to other solutions like the ones proposed 
in [14J among others, where other gravitational contributions might appear. The method 
could also be extended for the analysis of more complicated solutions like the Kerr one, 
discovered recently inside the massive gravity formulations ca For the case of solutions 
containing the time variable inside the metric components, the Vainshtein conditions (10) 
cannot ignore the time variation in such a case. Then the parameter combinations where 
the Vainshtein (screening) mechanism appears might change with time. 
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